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A THERMAL EQUATION FOR FLAME QUENCHING ' 

By A. E. Potter, Jr., and A. L. Berlad 


SUMMARY 

An approximate thermal equation was derived ff or quenching 
distance based on a previously proposed diffusional treatment. 
The quenching distance was expressed in terms of the thermal 
conductivity, the fuel mole fraction, the heat capacity, the rate 
of the rate-controlling chemical reaction, a constant that depends 
on the geometry of the quenching surface, and one empirical 
constant. 

The effect of pressure on quenching distance was shown to be 
inversely proportional to the pressure dependence off the flame 
reaction, with a small correction necessitated by the eflfect off 
pressure on flame temperature. 

The equation was used with the Semenov equation for burning 
velocity to show that the quenching distance was inversely propor- 
tional to burning velocity and pressure at any given initial 
temperature and equivalence ratio. 

In order to use the equation with experimental data, the rate- 
controlling reaction must be specified. Two choices were made: 
(1) the reaction between active particles and fuel ( corresponding 
to the difffusional equation), and (2) the reaction between oxygen 
and fuel. 

The two equations, one based on the active-particle— fuel 
reaction, the other based on the oxygen-fuel reaction, were tested 
using published data for the quenching off propane-oxygen- 
nitrogen flames. The data selected included the effect on 
quenching distance off oxygen-nitrogen ratio, propane concen- 
tration, and initial mixture temperature and pressure. A 
correlation off these data was obtained using each off the two 
equations, although both possessed shortcomings. The equation 
using the active-particle — ffuel reaction did not correlate data 
for all rich mixtures; the equation using the oxygen-ffuel reac- 
tion, while correlating data ffor both rich and lean mixtures, 
showed, a larger deviation ffrom the predicted linear relation 
than the other equation. 

INTRODUCTION 

The process of flame quenching is of interest, since it may 
be related to other combustion phenomena of engineering 
importance, such as flame, stabilization, flammability limits, 
and the general behavior of flames near cold walls. Flame- 
quenching processes become especially important in turbojet 
combustion systems when operation at low pressure (i. e., 


high altitudes) is considered. Flame quenching is usually 
studied experimentally in terms of the quenching distance, 
which is defined as the minimum channel size that will allow 
a given flame to propagate. Most of the available quenching 
data have been obtained either by (1) observation of the 
minimum tube diameter or rectangular slot width that will 
allow a flame to flash back, or by (2) determination of the 
minimum distance between plane-parallel plates that will 
allow a flame to propagate from a spark of minimum igni- 
tion energy. 

A rigorous theoretical treatment of the quenching process 
appears to be a most difficult task (refs. 1 and 2). Conse- 
quently, present theoretical treatments of quenching are 
necessarily approximate and seek primarily to correlate the 
quenching process with the variables that affect it, such as 
pressure, temperature, fuel type and concentration, inert- 
gas concentration, and quenching-surface geometry. 

Approximate treatments of quenching have been based 
either on a thermal or a diffusional mechanism for the process 
(refs. 3 and 4). Such apparently different outlooks yield 
results consistent with experimental data because the equa- 
tions for heat and mass transfer are formally identical, and 
because the thermodiffusivity and the molecular diffusivity 
are numerically nearly equal (ref. 5). Thus, each purely 
thermal approach to a combustion process has a diffusional 
analog. 

A quenching-distance equation that successfully correlates 
lean quenching data (refs. 4, 6, and 7) is the diffusional 
equation given in reference 4. Unfortunately,, the use of 
this equation in its present form is limited to stoichiometric 
or lean mixtures. The source of this difficulty appears to 
be the choice of the reaction kinetics. 

Essentially, the objective of this paper is to extend, if 
possible, the useful range of the quenching concepts of 
reference 4 to include hydrocarbon-rich mixtures. Because 
of the similarities between heat and mass flow, it is to be 
expected that a thermal analog equation can be derived that 
will correlate data equally as well as the diffusional equation 
of reference 4. It is felt that such an alternate thermal 
equation is more susceptible to changes in the reaction 
kinetics than the original diffusional model. The result of 
such changes is an extension of the useful range of fuels and 
oxidants that may be treated. 


Supersedes NACA TN 3398, "A Thermal Equation for Flame Quenching,” by A. E. Potter, Jr., and A. L. Berlad, 1955. 
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This report contains the derivation of a thermal quenching 
equation in which the rate-controlling reaction is not speci- 
fied. Two possible rate-controlling reactions are postulated, 
and the resulting two quenching equations are tested using 
published quenching data. It is shown that, if the oxygen- 
fuel reaction is assumed to be rate-controlling, quenching 
data are satisfactorily correlated for both rich and lean 
propane-oxygen-nitrogen flames. The pressure dependence 
of the quenching distance and the relation of quenching 
distance to burning velocity are also discussed. 

THEORY 

It is possible to obtain the thermal equation directly from 
the dift'usional equation (if the initial assumption is slightly 
altered), but it is felt that the derivation presented below is 
more understandable and allows the assumptions involved 
to be clearly seen. 

For a quenching equation, an initial postulate concerning 
the conditions under which a flame can exist inside a tube 
is necessary. A reasonable postulate is as follows: If the 
heat lost by the flame to the tube walls exceeds a critical 
value, the flame will be quenched. Since a definite amount 
of heat is produced in a given flame, an equivalent state- 
ment is that, in order for a flame to exist, the heat retained 
by it must be greater than a critical amount. It will be 
assumed that this critical amount is a constant fraction of 
the total heat produced in the flame. The criterion may 
be written as 


H t >FH t (1) 

(Symbols are defined in the appendix.) 

For the combustion of hydrocarbons, complete combustion 
of the fuel is assumed to occur for all mixture compositions. 
Consequently, the total heat produced by the process may 
be written as the product of the mole fraction of fuel in the 
unburned gas X r and the heat produced upon the disappear- 
ance of 1 mole of fuel A FI (as defined herein, AH is constant 
and equal to the heat of combustion to C0 2 and H 2 0 only 
for stoichiometric or lean mixtures; for rich mixtures, AH 
decreases as the oxygen concentration decreases). Thus, 
for hydrocarbon flames, equation (1) may be written as 

H r > FX f AH (2) 

The heat retained by the flame can be written in terms of 
heat capacity and average temperature of the flame in the tube 


The average temperature of the flame may be found if it 
is assumed that heat flow to the walls follows a differential 
equation of the form 


(T-T 0 )"=! 


(5) 


Here, q is the rate of temperature rise caused by chemical 
reaction, x is the thermodiffusivity and primes refer to 
differentiation with respect to x. 

The result of integrating and averaging equation (5) with 
the boundary condition T—T 0 at the wall is the equation 


T-T 0 = 


d 2 q 

G t x 


(6) 


In this equation d is a characteristic dimension of the tube, 
and G t is a constant associated with the geometric shape of 
the tube. Derivations of G t values for various geometries 
are given in reference 7. Equations (4) and (6) may be 
combined to give 


d 2 ^ FX f AHx 
G t ~ q 


(7) 


In order for a flame to pass down a tube, the tube size, 
given by d 2 jG t , must be equal to or greater than the right 
side of equation (7). Since the quenching distance is the 
characteristic tube dimension that just quenches a given 
flame, equation (7) can be written as a quenching-distance 
equation by removing the inequality sign; 


d 2 = 


FGiXjAHx 


(8) 


Then the thermodiffusivity x can be written as 

v «r RT 

PC v ,r 

The rate of temperature rise q can be written as 

RTAHw 

2_ ^pZy\T 


(9) 


( 10 ) 


Here, w is the rate of disappearance of fuel, molecules per 
cubic centimeter per second. 

Equations (8), (9), and (10) may be combined to give for 
the thermal quenching equation 


d 2 -- 


FGjNlirX f 

C v , T w 


( 11 ) 


H t =C v {T-T 0 ) (3) 

Equation (2) can then be written as 

C p ( T— T 0 ) > FX f AH 


The most important fact concerning this equation is that 
the form for w is not specified; any reaction may be chosen 
as rate-controlling. 

Equation (11) may be derived from the diffusional equa- 
tion if it is assumed that the thermodiffusivity and diffusion 


(4) 
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coefficient are equivalent, and if the criterion for flame prop- 
agation is assumed to be that the number of reaction events 
per cubic centimeter must be a critical fraction of the total 
number of events that normally occur in the unconfined 
flame. 

In terms of the active-particle — fuel reaction (the reaction 
specified as rate-controlling in the diffusional equation of 
ref. 4), the reaction rate w is given by 

w=k t a r '52c i ,T ( 12 ) 


If a collision-theory-type temperature dependence is as- 
signed to k t , equation (12) becomes 


w=B t a T ^Zc i ,r Tr Vz exp 

Combination of equations (10) and (13) gives 

d2 _ FG t NKrX f 


G v , r B t a r '^2,c i , r T r y exp 


(=* 

\RT, 


?) 


F , 

= B+‘ 


(13) 


(14) 


The choice of rate-controlling reaction is not limited to 
the active-particle — fuel reaction. For example, it is possible 
to follow Semenov (ref. 8) and assume the reaction in the 
combustion zone to be bimolecular and first-order with 
respect to fuel and oxygen. For this assumption, the real 
reacting species are obviously not fuel and oxygen molecules, 
but such an assumption is satisfactory if the concentrations 
of the reacting species are proportional to the fuel and oxygen 
concentrations. In this case, then, the rate of the reaction 
is given by 

w=k s a T b r (15) 


If a collision-theory -type temperature dependence is as- 
signed to k s , equation (15) becomes 


w=B s a T b r T r y exp 

The combination of equations (10) and (16) gives 

,, FGjNK r Xf F . 

C p . T B s a r b r T r .» exp (j^r) * 


(16) 


(17) 


In order to test equations (14) and (17), data for the 
quenching of propane-oxygen-nitrogen flames were used. 
These data points include the effect on quenching distance 
of oxygen-nitrogen ratio, fuel concentration and pressure 
(ref. 6) , and initial temperature (ref. 9) . 

In order to use equation (14), the assumptions made in 
reference 4 concerning the active-particle — fuel reaction were 
followed: T r was taken to be 0.7 T p , the partial pressure of 
fuel in the reaction zone was assumed to be one-half the 


partial pressure of fuel in the unburned gas, the partial 
pressures of the active particles (H atoms, OH radicals, and 
0 atoms) in the reaction zone were assumed to be 0.7 times 
their partial pressure in the adiabatic equilibrium flame, 
and E, was assumed to be 7 kilocalories per mole (ref. 10). 

In order to use equation (17), the partial pressures of fuel 
and oxygen in the reaction zone were taken as equal to 
their values in the unburned gas and E s was assumed to be 
38 kilocalories per mole (ref. 11). In keeping with previous 
usage (refs. 4 and 6), T r was chosen to be equal to 0.7 T F . 

The thermal conductivities were calculated for the un- 
burned gas mixture at 0.7 T P . The conductivities for 
oxygen, nitrogen, and propane were calculated at 1553° K, 
using the tables and data given in chapter 8 of reference 12. 
Conductivities at other temperatures were computed on 
the assumption of a %-power temperature dependence. The 
thermal conductivities of the gas mixtures were calculated 
as the sum of the mole fraction times the thermal conduc- 
tivity of each component, or, 


n='y~\XiK i (18) 

i 

It was felt that this simple linear mixing rule was adequate 
for systems largely composed of nitrogen and oxygen because 
of then- similar conductivities. 

The heat capacities were calculated for the unburned gas 
mixture at 0.7 T F by the same linear mixing rule used for 
the thermal conductivities. Individual heat capacities for 
oxygen, nitrogen, and pi’opane were calculated from data 
given in refei’ence 13. 

The equilibrium adiabatic flame temperatures and product 
compositions were calculated by the matrix method of 
reference 14 using the thermodynamic constants of reference 
14 and the heat of formation of propane given in reference 13. 

RESULTS AND DISCUSSION 
CORRELATION OF QUENCHING-DISTANCE DATA 

Inasmuch as the square of the quenching distance is 
predicted to be directly proportional to (eq. (14)) or to 
4* s (eq. (17)), values of \f/ t and \p a were calculated for the pres- 
sure, equilibrium adiabatic flame composition and tempera- 
ture, unburned gas composition and temperature, and tube 
geometry (assumed to be plane-parallel plates) which corre- 
spond to the quenching distances reported in references 6 and 
9 for propane-oxygen-nitrogen flames. Plots of d 2 against \p, 
and \p g are shown in figures 1 (a) and (b), respectively. The 
ranges of the variables included were: pressure, 0.1 to 1.0 
atmosphere; equivalence ratio, 0.33 to 1.90; mole fraction 
oxygen in the oxidant mixture, 0.21 to 0.70; and unburned 
gas temperature, 300° to 558° K. In figure 1 (a), it is shown 
that p, (which is entirely analogous to the diffusional equa- 
tion of ref. 4 except for the choice of quenching criterion) 
satisfactorily correlates the effect of these variables except 
for propane-rich mixtures where <p)> 1.2. 


L 
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Inspection of figure 1(b) reveals that 4> s (which includes 
the oxygen-fuel reaction as the rate-determining step, 
rather than the active-particle — fuel reaction) correlates not 
only the propane-lean, but also all the propane-rich quench- 
ing data. Thus, it appears that by proper choice of the 
rate-controlling reaction, the quenching concepts of refer- 
ence 4 may be extended in such a way as to include rich 
mixtures for the propane-oxygen-nitrogen system. 

The same reaction-rate constants were used for mixtures 
both rich and lean in propane. This may not be possible for 
all hydrocarbons; it may be necessary to use different rate 
constants for rich and lean mixtures, since it is generally 
believed that the reaction mechanism changes in going from 
a lean to a rich hydrocarbon flame. 

A mean line drawn through the data presented in figure 1 
reveals that the predicted linear relation between d 2 and i/y 
and \p s does not hold exactly. Rather, the data are best 
correlated if d 2 oc (4'i) OM oz{4' s ) OM . This difficulty is related, 
at least in part, to the choice of reaction mechanism, since 
the “best correlation” exponent is changed by a change in 
reaction mechanism. 

PRESSURE DEPENDENCE OP QUENCHING DISTANCE 

Examination of the equation for quenching distance 
(eq. (11)) reveals that all terms on the right side are inde- 
pendent of pressure except the reaction rate w. However, 
hydrocarbon flame temperatures change slightly with 
pressure because of changes in the extent of dissociation. 
Consequently, a pressure change indirectly affects tempera- 
ture- and composition-dependent terms in the equation. 
This effect is quite small, so that the pressure dependence 
of the quenching distance may be regarded as primarily 
dependent on the pressure dependence of the initial chemical 
processes in the combustion wave. 

The calculated pressure dependence of the quenching 
distance (including the indirect effect of pressure on flame 
temperature and composition) is compared in table I with 
the experimental values of reference 6. The average devia- 
tion of the calculated values from the observed values was 
found to be 10 percent for \p t and 17 percent for if/ s . Reference 
6 gives a value of 7 percent for the average deviation of 
values computed for the diffusional equation. 

Sy It is seen that i p s does not predict the pressure dependence 
as well as does ip t , except for the high-velocity flames at an 
oxygen fraction of 0.70. This is probably because the assump- 
tion of a simple second-order reaction between fuel and oxy- 
gen is a poor one. Concerning prediction of pressure depend- 
ence, a better choice of over-all order would be 1.7 rather 
than 2.0; this is in agreement with theoretical studies of the 
pressure dependence of burning velocity (ref. 12, p. 765), 
which indicate that the global reaction in most hydrocarbon 
flames ranges between first and second order. Such a choice 
would also tend to improve the linearity of the relation be- 
tween d 2 (observed) and f s (calculated). 

QUENCHING DISTANCE AND BURNING VELOCITY 


TABLE I.— COMPARISON OF OBSERVED PRESSURE 
DEPENDENCE OF QUENCHING DISTANCE 
WITH PREDICTED VALUES 


Oxygen 
fraction, a 

Equivalence 
ratio, <p 

Exponent n describing pressure dependence, 
d <x p-n 

Observed 

Predicted 

Diffu- 
sional * 

Thcr 

mal 

0.17 

0. 943 

0.90 

0. 88 

0. 83 

1.03 


1.000 

.89 

.85 

.76 

1.01 


1.340 

.84 

.75 

.75 

1.00 


1.530 

.71 

. 75 

.75 

1.00 

0. 21 

0. 738 

0. 85 

0. 86 

0. 80 

1.00 


.864 

.84 

.89 

.84 

1.03 


1.000 

.89 

.93 

.89 

1.03 


1.240 

.95 

.78 

.79 

1.02 


1.490 

.98 

.76 

.77 

1.00 

0. 30 

0. 566 

0. 76 

0.87 

0. 87 

1.01 


.662 

.93 

.89 

.88 

1.04 


.773 

1.06 

.92 

.90 

1.07 


1.000 

.98 

.94 

.90 

1.09 


1.380 

.93 

.86 

.84 

1.06 


1.903 

.74 

. 76 

.76 

1.01 

0.50 

0.476 

1.01 

0. 91 

0. 89 

1.06 


.544 

1.01 

.93 

.91 

1.07 


.680 

.96 

. 9G 

.92 

1.09 


1.000 

.93 

.97 

.92 

1.11 


1.358 

.91 

.95 

.90 

1.09 


1.940 

.88 

.82 

.81 

1.05 

0.70 

0. 345 

1. 12 

0.91 

0. 89 

1.06 


.395 

1.12 

.93 

.91 

1.07 


.494 

1.02 

. uo 

.92 

1.09 


1.000 

1.01 

.98 

.93 

1.09 


1.234 

1.07 

.98 

.92 

1. 11 


1.829 

.60 

.89 

.86 

1.08 


* Ref. 6. 


equation relating the two variables is derived on a purely 
diffusional basis. The thermal analog to this equation may 
be derived as follows : 

According to reference 16, the Semenov equation for the 
burning velocity may be written as 


2 k f W 

doPo^p, F 



(19) 


Multiplying equation 

(ID 

by equation 

(19), converting 

specific heats to molar heat capacities, and writing a„ and p 0 in 

the form 






NX f P 

RT 0 

(20a) 


Po- 

MP 

~RT„ 

(20b) 


yield 


d- 


RT, 

UP 


, 1 2K r K F F6 t Q^ (fle 

' V C v ,rC v , F W 


kRT \ m — 

dc p p) f w 


( 21 ) 


If it is assumed that W and w have the same pressure depend- 
ence, it follows from equation (21) that, at a given equiva- 
lence ratio, 


The relation of quenching distance to burning velocity has 
been discussed in references 3, 4, and 15. In reference 4, an 


U oc 


dP 


Oj-p-Cre+l) 


( 22 ) 



6 


REPORT 1264 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


From reference 6, it is seen that n decreases from about 
— 0.9 to —1.05 as the percent oxygen in the oxidant mixture 
is increased from 21 to 70. Simultaneously, the burning 
velocity increases from about 40 centimeters per second to 
the order of 300 centimeters per second. Consequently, 
equation (22) predicts that the exponent describing the pres- 
sure dependence of U should increase from about —0.1 to 

0.05 as U increases from 40 centimeters per second to about 
300 centimeters per second. This result may be compared 
with the experimental findings of Lewis (ref. 17) (which 
indicate an increase of about — 0.05 to 0.07 as burning veloc- 
ity changes from 40 to 300 cm/sec) to confirm equation (22) 
qualitatively. 

CONCLUDING REMARKS 

The fact that reasonable agreement exists between experi- 
ment and quenching concepts based on either heat or mass 
transfer indicates only that these processes have the same 
formal laws and that the thermal and molecular diffusivities 
are either equal or remain directly proportional for the data 
considered in this study. No evidence for the preponderance 
of either heat or mass transfer in quenching can be deduced 
from these results. The principal advantage gained by 
formulation of the quenching concepts of reference 4 in a 
“thermal” form is that a change in reaction kinetics is much 
simpler from a conceptual point of view for the “thermal” 
form than for the “diffusional” form. This conceptual 
advantage might be useful in a discussion of bizarre fuel- 
oxidant systems for example, CS 2 — F 2 flames. In addition, 
the actual calculations may be simplified, since by the proper 
choice of reaction mechanism, it may become unnecessary 
to calculate the equilibrium flame composition. 

SUMMARY OF RESULTS 

The results of an investigation into the possibility of ex- 
tending quenching concepts previously proposed to include 
quenching of hydrocarbon-rich flames may be summarized 
as follows: 

1. An equation for the critical quenching configuration of 
a channel was derived for hydrocarbon flames : 


^fgjNktXs 

C„, r W 

where 

C V , T heat capacity in reaction zone, cal/(°K) (mole) 
d characteristic dimension of tube geometry; quench- 

ing distance, cm 

F constant that relates total heat produced by com- 
bustion to heat which must be retained by flame 
for it to exist 

Gt dimensionless factor, dependent only on channel 
geometry 

N Avogadro’s number 

w rate of reaction in reaction zone, molecules/ (cc) (sec) 
X f mole fraction of fuel in unburned gas 
k t mean thermal conductivity in reaction zone, cal/ 

(cm) (sec) (°K) 

2. Two possible flame-initiating reactions were considered 
in the detailed formulation of w: (1) the reaction between 
active particles and fuel molecules, and (2) the reaction 
between oxygen and fuel molecules. 

3. The two equations resulting from the two reaction 
mechanism choices were tested using published data that 
included the effect of oxygen-nitrogen ratio, fuel concentra- 
tion, pressure, and unburned gas temperature on the quench- 
ing distance for propane-oxygen-nitrogen flames. These 
variables were correlated reasonably well for both rich and 
lean flames by the equation involving the oxygen-fuel 
reaction. The equation involving the active-particle — fuel 
reaction was satisfactory primarily for lean flames. 

4. It was concluded that, by the proper choice of reaction 
kinetics, the quenching concepts previously proposed may 
be cast in a form useful for the prediction of both propane- 
rich and -lean quenching data. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 6, 1954 


APPENDIX 


SYMBOLS 


The following symbols are used in this report: 
a fuel concentration, molecules/cc 

B Arrhenius constant 

b oxygen concentration, molecules/cc 

C p heat capacity, cal/(°K) (mole) 

C P , F average heat capacity, T 0 to T F , cal/(°K) (mole) 

Ci concentration of ith active particle, molecules/cc 
c p specific heat, cal/(°K) (g) 

Cp average specific heat, T 0 to T F , cal/(°K) (g) 

D diffusion coefficient, sq cm/sec 

d characteristic dimension of tube geometry ; quench- 

ing distance, cm 

E activation energy, cal/mole 

F constant that relates total heat produced by combus- 
tion to heat which must be retained by flame for 
it to exist 

6i dimensionless factor, dependent only on tube 

geometry 

AH heat released upon consumption of 1 mole of fuel 
by combustion process, cal/mole 
H r heat produced by chemical reaction in primary 
reaction zone, cal/mole 

II T total heat produced bv combustion of 1 mole of 
unburncd gas, cal/mole 
k rate constant, (cc) (molecules) /sec 

M average molecular weight of unburned gas, g/molc 
m molecularity of flame reaction 

N Avogadro’s number 

n exponent describing pressure dependence of quench- 

ing distance 

ni/u -2 moles of reactant per moles of product from stoichi- 
ometric equation 
F pressure, atm 

q rate of temperature rise, °K/sec 

R gas constant 

T temperature, °Iv 

U flame speed, cm/sec 

W average reaction rate in flame front as defined by 
Semenov, molecules/ (cc) (sec) 

w reaction rate in reaction zone, molecules/ (cc) (sec) 

X f mole fraction of fuel 

x distance, cm 

a mole fraction of oxidant in oxidant-inert mixture 

k thermal conductivity, cal/(cm)(sec)(°K) 

p density, g/cc 

tp equivalence ratio 

X thermodiffusivitv, sq cm/sec 

G 1 Nk t X ! 

& n . rr 1/2 —E, 

p , r^r^rT r GXp ^ 

GjNurXf 

' ‘ Cp,ra r ^c iiT T r 1/2 exp 

Subscripts: 

F flame 

i active-particle species 


o unburned gas 

r reaction zone 

s involves reaction of oxygen and fuel molecules 

t involves reaction of active particles and fuel 

molecules 
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